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A numerical study of the ghost Dyson-Schwinger equation in Coulomb gauge is performed and 
solutions for the ghost propagator found. As input, lattice results for the spatial gluon propagator 
are used. It is shown that in order to solve completely, the equation must be supplemented by a 
nonperturbative boundary condition (the value of the inverse ghost propagator dressing function at 
zero momentum) which determines if the solution is critical (zero value for the boundary condition) 
or subcritical (finite value). The various solutions exhibit a characteristic behavior where all curves 
follow the same (critical) solution when going from high to low momenta until 'forced' to freeze out 
in the infrared to the value of the boundary condition. The renormalization is shown to be largely 
independent of the boundary condition. The boundary condition and the pattern of the solutions 
can be interpreted in terms of the Gribov gauge-fixing ambiguity. The connection to the temporal 
gluon propagator and the infrared slavery picture of confinement is explored. 
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I. INTRODUCTION 



The ghost sector in nonperturbative studies of quantum chromodynamics (QCD) has invoked considerable interest 
over the last decade. One tool for studying nonperturbative QCD is the set of Dyson-Schwinger equations which, since 
one works in the continuum, is especially suited for discussing the infrared behavior of the Green's functions where 
Q_i' dynamical singularities may be present. Originally, in Landau gauge it was believed that the ghosts are unimportant 
to the Yang-Mills sector - the conjecture was that it is the three-gluon vertex that is responsible for confinement 
[l|, 0|- However, the work of Ref. Q turned this conjecture on its head: it was shown that the ghost sector dominates 
the set of Dyson-Schwinger equations in the infrared, with the three-gluon vertex contributions being subleading. 

In order to obtain numerical results, the authors of Ref. [3[ introduced infrared fit functions for the propagator 
^ ' dressing functions that were characterized by powerlaws (in other words, ~ x K for some exponent k). Drawing 
on this idea, it was subsequently shown that one can completely characterize the infrared behavior of Yang-Mills 
theory in Landau gauge by studying only the ghost sector contributions, with a tree-level ghost-gluon vertex and 
pure powerlaws 0, |5[ . The powerlaw theme was expanded [6 8] and it was later shown that the tree- level ghost- 
• ■ gluon vertex truncation was reliable 9]. Two contemporary reviews can be found in Refs. [id lll|. These results 
gained some degree of popularity not simply because of their simplicity, but also because they were in agreement with 
various descriptions of the confinement problem: the results imply positivity violation, such that the propagators 
cannot correspond to physical states and formalized by the Oehme-Zimmermann superconvergence relations [12J, the 
results were also in agreement with the Kugo-Ojima [13[ and Gribov- Zwanziger [Tll - fl6j confinement scenarios (at 
least insofar as they were understood at that time). The lattice results of the day were in good agreement with the 
powerlaw description of the infrared (for a detailed discussion, see for example Ref. [ll]). However subsequent lattice 
studies going further into the infrared cast doubt on and eventually seemed to rule out altogether [lj, LL8| the existence 
of infrared powerlaw solutions for the propagator dressing functions in Landau gauge and indicated that the functions 
should be a finite constant at zero momentum (for a detailed discussion, see Ref. [19j ). In addition, it was pointed 
out that non-powerlaw solutions also exist to the Dyson-Schwinger equations (see, for example, Refs. [20, [21] and 
discussions therein). The coexistence of both types of solution to the Dyson-Schwinger equations in Landau gauge 
was studied in Ref. (22|. The issue at stake is how the confinement mechanism manifests itself at the level of the 
propagators in Landau gauge and the investigation is ongoing. 

The situation in Coulomb gauge studies of nonperturbative QCD is in some respects similar. Coulomb gauge 
studies come in various guises and indeed, the definition of Coulomb gauge is somewhat different depending on which 
formalism one uses. Currently, the most widely used continuum formalism is the canonical formalism, based on the 
Hamilton densit y op erator [23rl29| . In the canonical approach, one starts by setting Weyl gauge and then 'rotates' into 
Coulomb gauge [30], resolving Gauss' law to ensure gauge invariance and subsequently minimizing the energy density 
with an Ansatz for the wavefunctional to arrive at a set of Dyson-Schwinger-like equations. In the canonical formalism, 
there are no temporal degrees of freedom and the equations refer only to the equaltime correlation functions. The 
Dyson-Schwinger-like equations closely resemble the Landau gauge Dyson-Schwinger equations in a space with one 
less dimension. In the canonical approach, there is known to be two types of solution for the propagator dressing 
functions [27] (just like in Landau gauge): critical and subcritical and in the critical case, two different values for the 



powerlaw exponent have been found. The favored solution is the most singular solution where the ghost propagator 
dressing function diverges as l/|fc| in the infrared (\k\ — > 0) because this solution has a natural interpretation with 
regards to physical confinement (see the discussions of Refs. [28|, [29|). The corresponding expression for the dressing 
function of the equaltime spatial gluon propagator is very close to Gribov's famous formula [l4| (see later in the text 
for the explicit form). 

There also exists lattice results for Coulomb gauge [31-35]. Again, the definition of Coulomb gauge plays a role. 
In Refs. [31H33J, the gauge is implemented separately on each time-slice and the equaltime correlation functions are 
considered (as in the canonical formalism). However, these results appear to be plagued by scaling violations. Such 
scaling violations highlight one of the difficulties inherent to Coulomb gauge, namely that as yet there is no complete 
proof of renormalizability (although there has been some progress in this direction [l6|, [3a, [37|). In an attempt to 
avoid the scaling violations, the authors of Refs. 34, 35] implemented a lattice version of Coulomb gauge that includes 
the temporal degrees of freedom, thereby studying the full correlation functions and not their equaltime counterparts. 
All the lattice studies so far agree on certain results: i) the equaltime spatial gluon propagator is well-reproduced by 
the Gribov formula, ii) the ghost propagator is infrared enhanced with the dressing function diverging like l/[fc 2 ] K and 
where k sw 0.2 — 0.25 (at least for as far as these lattice studies go into the infrared regime), iii) the temporal gluon 
propagator is strongly infrared enhanced, probably diverging like 1/k 4 . We should however remind the reader of the 
earlier situation for Landau gauge studies (discussed in Ref. [19]) - there was a generally agreed infrared result for 
modest lattice sizes (sizes comparable with the current Coulomb gauge lattice studies) which was later overturned by 
much larger lattices that probed further into the infrared region. Therefore, one should perhaps take some caution in 
interpreting the lattice results in isolation. However, one point of note is that the Gribov formula for the equaltime 
spatial gluon propagator in Coulomb gauge vanishes in the infrared and this was shown to be true in general on the 
lattice [38} (though there are certain caveats with this, see for example Ref. J3l} for a discussion). Taken alongside 
the results of Ref. [38| , of the continuum canonical approach [2J, [26| (discussed previously) and in the absence of any 
direct evidence to the contrary it seems fair to say that there is at least a current agreement about the equaltime 
spatial gluon propagator having the Gribov form (whether this situation will persist may well prove interesting). 
Certainly, it appears to us reasonable to use the Gribov formula as input into the work presented here. 

A third approach to Coulomb gauge (and that used in this study) is the continuum functional formalism which is 
based on the Lagrange density of QCD. The importance of the functional formalism in Coulomb gauge lies in the 
recognition that one can reduce the system to physical degrees of freedom [l6[. Further, one can formally show the 
existence of a conserved and vanishing total color charge (along with the absence of the infamous Coulomb gauge 
energy divergences) [39j. This is crucial because a system that spuriously 'leaks' color charge cannot be confining. 
Because of the inherent noncovariance of Coulomb gauge and the fact that one explicitly retains the temporal degrees 
of freedom in the functional formalism, detailed technical results are difficult to derive. However, steady progress 
is being made: the Dyson-Schwinger equations have been explicitly derived [40h42| . along with the Slavnov- Taylor 
identities [4l|, [43|, [44j and one-loop perturbative results are available [4l|, [42|, |45| . These results include a study of heavy 
quarks and the corresponding Bethe-Salpeter equation [44j that motivates an extremely simple quark confinement 

scenario that corresponds to the infrared slavery picture whereby if the temporal gluon propagator diverges like 1/k 4 
in the infrared (as indicated from lattice results, discussed previously) then one has only colorless finite energy bound 
states of quarks and antiquarks with a linearly rising potential between them. 

In this study, we shall investigate numerical solutions to the Coulomb gauge ghost Dyson-Schwinger equation 
within the (second order) functional formalism and focus on the issue of critical versus subcritical infrared behavior. 
As primary input, we utilize the general consensus about the equaltime spatial gluon propagator having the Gribov 
form, as discussed previously. The paper is organized as follows. In Sec. HH the ghost Dyson-Schwinger equation and 
its components will be introduced. In Sec. IIII1 the reduction of the equation to a form suitable for numerical analysis 
will be made explicit, along with a discussion of the asymptotic behavior. The numerical results will be presented in 
Sec. IIV1 the comparison of the solutions with available lattice results will also be described. In Sec. [V] the results will 
be discussed and arguments put forward about their connection to the Gribov gauge-fixing ambiguity. In Sec lVIl it 
will be motivated how the critical ghost solution gives rise to a straightforward interpretation of confinement. The 
paper closes with a summary and outlook in Sec. IVIII 

II. THE GHOST EQUATION 

Let us begin by reviewing some basic results for the Green's functions in the standard second order functional 
formalism as applied to Coulomb gauge Yang-Mills theory. Throughout this work, we shall use the notations and 
conventions established in [4fj, |4l|, |43j. We work in Minkowski space with metric g^ — diag(l, —1) (until such time 
as it is necessary to analytically continue to Euclidean space). Roman subscripts (i,j, . . .) denote spatial indices (all 







FIG. 1: The unrenormalized Dyson-Schwinger equation for the ghost two-point proper function, T- c . Filled blobs denote 
dressed two-point functions and empty circles denote dressed proper vertex functions. Wavy lines denote proper functions, 
springs denote connected (propagator) functions and dashed lines denote the ghost propagator. 

minus signs associated with covariant/contravariant vectors are explicitly extracted) and superscripts (a, b, . . .) denote 
color indices in the adjoint representation (with N c colors). 

The unrenormalized ghost Dyson-Schwinger equation is given by [4l[ (see also Fig. [TJ) 

T<£{k) = ? ^fc 2 + J duT^tik^- k, -«)W£(* - oj)WT Aij (u;)Tf c e Aj (k - w, -k,u), (2.1) 

where dw = d d+1 ui/(2ir) d+1 and with spatial dimension d. The ghost propagator (Wz c ) and the ghost proper 
two-point function (Tec) are decomposed as follows: 

Wg(k) = -S ab ^-D- cc (k 2 ), Tt c {k) = S ab ik 2 T- cc (k 2 ), (2.2) 

k 2 

and the respective (dimensionless) dressing functions obey the usual relationship: 

D- cc {k 2 )T- cc {k 2 ) = 1. (2.3) 

Notice that in general, Coulomb gauge propagators and vertices are not dependent on the (Lorentz invariant) four- 
momentum squared (fc 2 = k 2 . — k 2 ), but rather on the energy squared (feg) and the momentum squared (k 2 ) separately 
because of the inherent noncovariance of Coulomb gauge. However, we do know that the ghost propagator dressing 
function, D^c, is strictly independent of the energy as a nonperturbative result of the Slavnov- Taylor identities [4lL l43jj . 
In addition to the ghost propagator, we have the spatial component of the gluon propagator 

Wf AlJ {k) = 6 ab ^t l3 (k)D AA (klk 2 ), (2.4) 

where Uj(k) — 8ij — kikj/k 2 is the transverse spatial projector. The above dressing functions (Z^o D AA ) reduce to 
unity at tree- level and the explicit one-loop perturbative expressions are known [4l|, |45[ . The third component of the 
Dyson-Schwinger equation, Eq. (|2.ip . is the spatial ghost-gluon vertex which is decomposed as follows 

IS&(P1,P2,P8) = -Wf abC TccAi(pi,P2,P 3 ), (2.5) 

and at tree-level 

r HoL(^i'P2,P3) =Pu- (2.6) 

The spatial ghost-gluon vertex has the property that it reduces to the bare vertex in the limit of vanishing 'in-ghost' 
spatial momentum (\p 2 \ — » 0) |4Q] (as in Landau gauge [46]). With these decompositions, the ghost Dyson-Schwinger 
equation, Eq. (|2.ip . can be rewritten in terms of the dressing functions (still unrenormalized and in Minkowski space): 

r (t2\ i 2 Ar f duiD AA (oj 2 ,uJ 2 )Dc C ((k-uj) 2 ) 

Tcdk ) = 1 - ig N c / — — hUj{wWccAj{k ~ w, -k, w). (2.7) 

J k z uj z {k — u>Y 

We immediately notice that in order that the loop integral be independent of the external energy (&o), the spatial 
ghost-gluon vertex (or more properly, its contraction within the integral) can only be dependent on the energy of the 
gluon leg (u> ), i.e., 

TccAj{k~ oj,—k,u) -»■ Tcc A j(k — u,-k;u>). (2.8) 



This is entirely consistent with the vertex Slavnov- Taylor identities in Coulomb gauge [43J. Further, we notice that 
the infrared limit of the Dyson-Schwinger equation coincides with the known infrared limit of the dressed vertex. This 
allows us to specify the truncation scheme whereby the fully dressed vertex is replaced with its tree-level counterpart: 

FccAj(k-3, —k',u) -> (k — u)j. (2.9) 

This truncation is known to be robust in Landau gauge studies [a, lZJ, since it reproduces both the infrared and 
ultraviolet behavior of the vertex; the dynamical corrections to the dressed vertex in the mid-momentum region have 
been demonstrated to be very modest Q • Because the ghost couples in both Landau and Coulomb gauge to transverse 
gluon degrees of freedom, we can reasonably assume (though this should of course be verified at some stage) that this 
truncation will work well in Coulomb gauge. Indeed, it has been successfully employed in the analyses performed in 
the canonical approach [24], [2g, [27| ■ 

Before continuing, let us briefly introduce some aspects of the renormalization that are relevant specifically for the 
ghost Dyson-Schwinger equation. We assume that the functional approach to Coulomb gauge is (nonperturbatively) 
multiplicatively renormalizable in this study, although as mentioned in the introduction, there is currently no complete 
proof. When renormalizing the theory, the process of regularization introduces a nontrivial scale. Recalling the 
noncovariant nature of Coulomb gauge, we initially choose to renormalizc at the purely spatial momentum point 
feg = 0, k 2 = [i (where /i is finite) and assign the fields A, c, c (we will discuss the temporal component of the gluon 
field, a = Aq, later) and coupling (g) the renormalization coefficients y ' Za, \[Z~ c , and Z g (the ghost and antighost 
fields share the same coefficient). We define the renormalized coupling via the spatial ghost-gluon vertex; that means 
that in the effective action term for this vertex ~ gccA with the corresponding renormalization coefficient Za, we 
identify Za = Z g Z c \J Za, or 

Z g = — %=, g = gZ g (n;g), (2.10) 



where g is the renormalized coupling (and is a function of /i in the sense that with a different renormalization point, 
one would extract a different value for the coupling) and the renormalization coefficients are all functions of [i and 
~g. We leave the question of the regularization procedure (and possible scale) to one side for now. For the propagator 
dressing functions we then have 

D A A(k 2 ,k 2 ;g) = Z A (fi;g)D^ A (k 2 ,k 2 ; l i,g), D 5c = Z c Dg (2.11) 

(we shall leave aside the common arguments for notational convenience where appropriate). Conventionally, all 
renormalized propagator dressing functions are defined via the renormalization point: D R {k^ — 0, k 2 — /i; n,g) = 1 
(although this will not be explicitly required in this work). As mentioned above, we define the renormalized coupling 
via the spatial ghost-gluon vertex. To be more precise, we demand that at the renormalization point (p® = 0, P2 = 0, 
Pi = P% = M) the renormalized vertex is bare, i.e., 

r|U(p§ = o,P2 = o,pI =pI = mii,g)=pu. (2.12) 

However, the unrenormalized vertex is also bare at this point and we have that 

Z A (^g) = l =>Z g y/zZz e =l. (2.13) 

Clearly, this is the same situation as in Landau gauge pa ]. 

As mentioned in the introduction, the spatial gluon propagator in Coulomb gauge takes definite meaning only when 
one specifies within which formalism one is working. In the second order functional formalism used here, because 
there is no gauge restriction on the temporal component of the gluon field, the spatial gluon propagator retains 
its full temporal dependence. The Fourier transform associated with the complete spatial gluon propagator reads 
(dk =d d k/(2n) d ): 



WAAij(x°,xi]x^,x 2 ) = / ake- %k °(*°-*% +ti -&- Sa )W& ij (k) 



S a b f d k e M^-^) t .. { % )t r ^o e -»fc (*;-*°) D AA {klk 2 ) 



UAk)i / Zie-*^-"!) ^-^ — '—. (2.14) 

3 ' Loo 2tt (k 2 -k 2 + iO + ) V ' 

In the above, we have made use of the translational invariance and inserted the Feynman prescription for dealing with 
the poles of the propagator. Note that the multiplicative renormalizability of the theory (discussed above) refers to 



the complete propagator. Now let us further consider the equaltime case x® = x§. This is the spatial gluon propagator 



that one considers within the canonical (Hamiltonian) formalism 2J] and on the lattice where Coulomb gauge is fixed 
separately on each time-slice |3ll — [33| . Thus, we require that the integral (which will play a crucial role in the ghost 
self-energy studied in detail in the next Section) 

d t a{ p )=1 r dk ° DAAjkiP) = r° dk A D AA {-kiP) (2i5) 



„ 2vr ( fc 2 _ k 2 + l0+) J_ oo 2n (fc| + fc2 ) 

(the second form arising from the Wick rotation: fco — > tk^, which we assume to be possible) be well-defined. This 
places a weak restriction on the large (Euclidean) energy behavior of the dressing function D aa for arbitrary, finite, 
momentum k 2 . One important observation about this restriction concerns the perturbative limit and highlights a 
difficult technical problem associated with the functional approach to Coulomb gauge. Perturbation theory is only 
valid in the vicinity of the renormalization point (chosen here to be purely spacelike, i.e., k$ — 0); away from 
this point, the corrections grow logarithmically. When integrating over the energy, one is extending into the high 
(Euclidean) energy region where the perturbative result is invalid. Thus for example, in order to compare a numerical 
result for Daa with the (large momentum) perturbative expression for D T AA in the canonical approach, one already 
requires nonperturbative information about Daa, although the asymptotic series expansions still agree at one- loop 
[47J . Similar observations apply for the temporal component of the gluon propagator [48| . 

In the next section, it will be seen that in order to solve the ghost Dyson-Schwinger equation, one requires D T AA 
as input, as opposed to Daa- As was discussed in the introduction, there is a current consensus that D AA takes 
a form consistent with the Gribov formula. For definiteness however, let us begin by discussing the results of one 
particular lattice study of Daa [34] (to our knowledge the only calculation of this quantity to date) . These results 
can be summarized as follows: 



D AA (-ki,k 2 ) = 



1 



K 



i a (gi) 



VI 4 



(2.16) 



The lattice data are in Euclidean space, with the gauge group SU(2) and with a lattice coupling gi. The scale 
m sw 0.88GeV ~ 2^/5^ where a w is the Wilsonian string tension. It is observed that m is not dependent on the 
renormalization scale. The exponent a(gi) has the following behavior: 

, s f 0, qi 'small' ,_ _, _ N 

a ^) ™{i% 'large' " ^ 

Importantly, the temporal behavior of the dressing function factorizes into a dimensionless function and one obtains 
with Eq. (|JT5|) 

T p F(l/2-a(,g,)) 1 VP 

DAAk ] ~ r(i/2)r(i-«( 9; ))27f^T- (2 - 18) 

The temporal behavior of Daa collapses to an overall constant prefactor which will be seen in the next section to 
be largely irrelevant. The restriction, following the definition of D AA , Eq. (J2.15I) . on the large (Euclidean) energy 
behavior of Daa is related to whether or not this constant prefactor is finite or divergent. Clearly, the 'small' coupling 
result: a(gi) — is more pertinent since we are interested in the physical regime where the (renormalized) coupling 
is small and for 'large' coupling, the system undergoes various nontrivial phase transitions on the lattice. Let us thus 
write 

D T AA (k 2 ) = \ J^— (2.19) 

2 V fc 4 + m 4 

such that D AA — > l/(2|fc|) as k — > oo as an arbitrary normalization condition for now (it corresponds to the tree-level 
result in the canonical formalism [24], |47| ) ■ 

The expression Eq. (|2.19p is of course Gribov's original formula for the equaltime spatial gluon propagator [14[ 
and as discussed in the introduction, there is a current consensus that this is the correct (or close to correct) result. 
Notice that the above lattice result is already renormalized and there is no reference to the renormalization scale 
(as mentioned, the Gribov scale m is observed to be independent of the renormalization scale) except implicitly 
through the normalization. Interestingly (and unlike in Landau gauge), for large momenta there is no evidence on the 
lattice for a perturbati ve g luon anomalous dimension within D AA , i.e., the coefficient of the perturbative logarithms 
is consistent with zero 34] and this is consistent with the other lattice studies I3ll-l33ll. 



III. ANALYTIC FRAMEWORK 



Inserting the truncated form for the dressed ghost-gluon vertex, Eq. (J2.9I) , and rewriting in terms of renormalized 
two-point dressing functions, the ghost Dyson-Schwinger equation, Eq. (|2.7[) now reads 



Tg(k 2 ) = Z c -g 2 N c 



duJDg{{k-ujf) 



-/,',,:-, y -1 v I 'c. u/x^.v.yn.-uy J j k . ,r,\ f°° duJ D% A (u>Q,ui 2 ) 

p(fc-w)2 4 3 tA ' J-oo ^ H~^ + iQ + y 



(3.1) 



Because the truncated form of the vertex is not dependent on the energy, we immediately recognize D AA , given by 
Eq. (|2.15p and inserting our lattice input, Eq. (|2.19p . we then have 



rf c (fc 2 ) 



:9 2 N C 



du}VuPD£((k-uJ) 2 ) 

k 2 Vtu 4: + m 4: (k — Cj) 2 



K{Kjtij {bJ)- 



(3.2) 



Recall that D AA was arbitrarily normalized to remove the coefficient (the combination of gamma-functions) resulting 
from the energy dependence of D AA . Were we to relax this normalization condition, we would have some overall 
(constant) prefactor for the ghost self-energy. 

Before continuing, let us consider the perturbative treatment of Eq. (|3.2[) at one-loop (obtained by setting m = 
and Dec within the integral to unity). Using dimensional regularization, we get the standard result [4l| 



R(ti\ 



r£(* 



A g 2 k 2 

3 {An) 2 \ n 



A a 2 

-N — 

3 c (4tt) 2 



1 7 

7-7+ o 



21n2 + ln(47r) 



(3.3) 



from which one can identify (A = g /(Air) — a s /(Aiv)) 



Tg(k 2 ) = l + ^N c ln(j)+0(g i ) 



(3.4) 



At this level in perturbation theory, one need not consider the running of the coupling (i.e., g = g). Assuming that 
a s is small and k 2 is close to /i so that the logarithm is also small, the leading order expression can be resummed and 
inverted to give 



7s 



DUk 2 -^,g) 



/'■ 



7s 



-XN r 



(3.5) 



where 7 9 is the leading order expression for the ghost anomalous dimension. 

As it stands, Eq. (I3.2[) contains logarithmically UV-divergent pieces. To eliminate the divergences, we use a 
nonperturbative subtraction. Taking the finite renormalization scale [i as our subtraction point, we then have 



R/f2 



r£(Jr) = rg(lr = „) 



\fN, 



duj VlS 2 I Dg((k - UJ) 2 ) kikj 



V^ 4 + 



(k-ti) 



7fi2 



k 2 



tij{S) _pm^i^ Uj{S) 



(k-uj) 



:r\2 



k 2 



(3.6) 



fc 2 =i. 



One can see that (as long as D? c (x) is smoothly varying for k 2 , /i <C x — >• oo) the integrand then falls faster than l/|w| 3 
at large |cD| to ensure convergence. Translating the integration variable lu — > k — uj, the integral can be conveniently 
rewritten using a UV-cutoff (A) and introducing some notation, we write: 



x 

tfeu 



y = w , 



k-uj 



xyz, 



(An)' 



/ dz dy^/y, A -> oo, 



(3.7) 



to give (we reinsert the renormalization scale dependence of the functions and drop the sub- and superscript notation 
for T~ c for clarity) 



r(ar; fx) = T(^; fx) - XN C / — T(y; ^y 1 [I(x, y; m) - IQj,, y; m) 
Jo V 



(3.8) 



where the (angular integral) function / reads 



I(x,y;m) 



dz(l - z 2 



1 



x 

I) 



2j X -z 



-1/2 



(1 



V 



y 



-1/2 



(3.9) 



Note that after the UV- divergence has been subtracted, the above representation of the self-energy integral is exact, 
as long as A — > oo. Also, when the Gribov scale m = 0, one recovers the free gluon propagator. For m — 0, the 
angular integral can be performed analytically and 



I(x,y;m = 0) = - 







3/2 



(*-»)(-) +e(y-x) 



(3.10) 



With this, one recovers the perturbative result, Eq. (|3.4[) . when expanding in a s : 

T(x;n)=r(n;n) + ^\N c hi(-\ 



(3.11) 

At this stage, one can make a preliminary infrared analysis of Eq. (|3.8[) . For fixed g 2 and 771, given that the dressing 
functions are dimensionless, one can make the initial Ansatz that T(x; fj,) = T(x/n), disregarding for now the influence 
of the Gribov scale m. One can also infer from the sign of the perturbative logarithm that T(x//j,) increases with x, 
discounting for now the possibility that the function has some turning point (this will be seen numerically not to be 
the case). Thus, we can write for x/fi —> 



T(x/fi) = 00 + Ci - 



flO: a. 1 > 0, a > 0. 



(3.12) 



The constant ao is determined by the constant part of Eq. (|3.8|) . To get the powerlaw solution, we assume that this 
constant vanishes (meaning that T(x/ fi) ~ (x/ /i) a ) and rewrite Eq. p.8[) in terms of x/fj, and y / ' /i: 

T(x/n) = T(l) - XN C f ^Lviy/n)- 1 [I(x/n, vim m/sfH) - 1(1, y/n; m/VM)] - (3-13) 

Jo V 

noting that the function / is dimensionless. The integrals (angular and radial) will not alter the exponent of the infrared 
expansion since they are independent of x and /^, so to get the leading powerlaw relationship for the exponent, one 
must simply remove dy/y and replace y//i in the integral with x//x, similarly removing dz and setting z — (in other 
words, dimensional analysis!) and then let x/ jjl —¥ 0. Explicitly 



T{x/fi) ~ T{x/n) 1 [I(x/n, x/n; my/JJ) - 1(1, x/n; m^/JL)} 



where 



I(x/n, x/fi', m/y/JI) - I(l,x/fx; rriy/jl) 
so that the powerlaw relationship reads 



3/2' 



-1/2 



-Q + l 



(3.14) 



(3.15) 



(3.16) 



or a = 1/2 if one counts the powers of x (— k 2 ). If one counts instead the powers of fi, then one has a = which is 
the constant solution. To summarize, we would expect that if the dressing function T vanishes in the infrared, there is 
a leading powerlaw behavior given with the exponent a = 1/2. This is in agreement with the results obtained within 
the canonical formalism [2J, [2f| [27| ■ We will see later that the numerical result (derived independently of this crude 
analysis) confirms this behavior. 

Returning to Eq. (|3.8j) . for m / the ^-subtracted equation is still not suitable for numerical analysis because for 
general r~ 1 (y;/i) within the integral, the integral can become negative - this is obvious from the (not resummed) 
perturbative result, Eq. (|3.4p . (for which r~ 1 (y; /i) = 1) when x <^ /i, where In (x/ fx) — > —00 . Therefore, an iterative 



1 Demanding that such a negative integral be explicitly excluded perturbatively led to the original derivation of the Gribov factor used 
here as input tl4l| . 



8 



procedure to solve will automatically fail since we require T~ 1 (y;ii) to be finite and positive within the integration 
domain (it can have at most an integrable singularity at y — 0). To avoid such a situation, wc subtract Eq. (|3.8[) once 
more but now at the position x = to get 



r(ar; n) = r(0; /*) + XN C 



dy 

V 



r (2/;M) 1 [I(0,V,m) -I(x,y;m)] 



(3.17) 



Note that the above equation is still formally renormalized at the scale /i ^ 0. The bracketed combination of angular 
integrals is however positive definite for all x > (and zero for x — 0) . The dependence of the above equation on the 
original renormalization point is now only implicitly given within the coupling (and the normalization condition for 



D AA ). Further, one can rescale the variables x 



xra 



y — > ym such that the new (dimensionless) x, y are in units 



of [GeV] 2 /m 2 . This is a consequence of the fact that the form of the renormalized gluon propagator was fixed from 
the beginning. The rescaling proportional to m 2 also means that the comparison with the leading order perturbative 
result (where m — 0) will at most be for asymptotically large x. Thus, for fixed a s — 4n\ as input, the solution is 
independent of jjl and we can write 



r(a;) = T(0) + XN C 



dy 



w 



V 



3^1 



I(x,y;l) 



(3.18) 



where we have evaluated 1(0, y; 1). There are now only two input parameters, AiV c and T(0). For definiteness, we 
choose the physical values a s (M z ) — 0.1187 and N c — 3 and consider the solution with a range of values for T(0). 
From the above equation, it is now obvious that the overall normalization of D AA can be absorbed into an irrelevant 
constant prefactor for r(x). (This observation also applies to the comparison of the case N c — 2 with N c — 3.) Thus 
we observe that the only relevant information about the original input gluon propagator is that it had the Gribov 
form at equal times, the values of the Gribov scale and normalization condition will not affect the overall features of 
the ghost solution. The numerical solutions to Eq. (|3.18[) will be discussed in the next section. 

Having discussed the renormalized ghost Dyson-Schwinger equation, let us now briefly discuss the associated renor- 
malization coefficient Z c (p,,~g), Inserting the form for the spatial gluon propagator dressing function, Eq. (|3.2|) can be 
written in terms of the angular integral and the rescaled variables in the following way 



Z c (A,[a s ,r(0)]) = r(x) + AiV c 



dy 

yT(y) 



I(x,y,l). 



(3.19) 



In the above, we have recognized that Z c is dependent on A, A = Ai:a s and T(0). The renormalization scale (p.) 
dependence is implicit within a s , which is here fixed. If the Dyson-Schwinger equation is multiplicatively renormal- 
izable, then Z c should be observed to be independent of x and this is numerically verified (the deviations are of the 
same order as the tolerance for the convergence of the iterative solution to Eq. (J3.18I0 . It is also seen that Z c is 
independent of the numerical IR cutoff. Since Z c is independent of x, then we could equivalently set x = and use 
the simple form 



Z c (A,[a s ,T(0)})=T(0) + UN c [ dy 

3 Jo ^l + y z T(y) 



(3.20) 



IV. NUMERICAL RESULTS 



The results (for T(x) and D(x) = T(x)~ 1 ) are shown in Fig. [2] along with a fit to the supposed infrared powerlaw 
and the asymptotic perturbative result. Numerically, the solution is stable apart for the deep infrared region in the 
case of small T(0). In this case one has an integrable singularity and we explicitly do not include the deeper level of 
numerical sophistication required to deal with the problem since this would necessarily involve an assumption about 
the infrared behavior of the solution. However, the features of the solution are clearly evident. 

In the UV region, all solutions with different T(0) < 0.1 tend to the same curve (the case T(0) = 1 will be discussed 
presently). Asymptotically, it is seen that these solutions obey 



r(x)- l = D£ c (x) = 



p7g 



(4.1) 



where c « 1.581 and confirming that one recovers the (resummed) perturbative behavior for large x. That all curves 
have the same asymptotic large x behavior suggests that for a large renormalization scale fi (and after scaling variables 
with m 2 , n would have to be asymptotically large), there is a unique value of T(fi) (or vanishingly small differences 
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FIG. 2: Results for the ghost dressing functions F(x) [left panel] and D(x) [right panel]. See text for details. 
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A 


10 6 


10 8 


10 6 
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10 6 


10 8 


10 6 


10 7 


10 8 


10 9 


z c 


1.45 


1.57 


1.09 


1.24 


1.09 


1.24 


1.09 


1.17 


1.24 


1.31 



TABLE I: Numerical results for the renormalization coefficient Z c . Dimensionful quantities are in units of [GeV] 2 /m 2 , see text 
for details. 



between different values) for arbitrary input T(0). We can thus infer that our input T(0) (at least for moderate values) 
has nothing to do with the perturbative regime and its renormalization - it is genuinely nonperturbative in origin. 

Continuing with the UV region, results for Z c are given in Table|TJ Effectively, Z c ~ 1 and for the A and r(0) values 
sampled is slowly varying, increasing with A like Z C (K) ~ In A as one would have expected from the perturbative 
result (and with fixed spatial gluon input and a tree- level vertex, this is not surprising). Noticing the behavior of Z c 
for different values of input T(0), we plot Z c as a function of T(0) for fixed A = 10 6 [GeV] 2 /m 2 in Fig. [3] It is seen 
that Z c is roughly constant for input values less than T(0) ~ 0.3 and this gives an estimate of the largest value of 
r(0) for which the ghost solution does not change in the UV. This maximum value of T(0) loosely corresponds to 
the situation where the solution curve T(x) loses contact with the common (asymptotic) large x curve in Fig. [5] (as 
illustrated for the case T(0) = 1). The inference is that there is a region of values < T(0) < 0.3 for which the input 
constant T(0) is unimportant to the UV. We shall discuss this later in more detail. 

Returning to Fig. [51 in the infrared region for small T(0) (< 0.1) the solution has a region for which the naive 
powerlaw solution holds and as T(0) is decreased, one can see that this region extends further into the infrared. 
Moreover, extrapolating to the case T(0) = by eye, one can see that the pure powerlaw solution does exist. Let 
us emphasize that the powerlaw has not been artificially introduced into the numerical solution via sophisticated 
techniques (and for which the numerical stability could be improved to an arbitrary degree) . The results are obtained 
using iteration and a direct numerical integration grid: the observed infrared behavior is an independent result. The 
'fitted' (i.e., fitted by eye) form of the powerlaw solution has one input parameter - the coefficient (the exponent was 
fixed from the earlier analysis). Clearly, the coefficient must have the factor \fXW c but that the factor y/2 works so 
well is not explained (and given the previous discussions about the normalization of D^ A is not really relevant). To 
summarize the IR behavior: the most notable feature of the various solutions is that when going from the common 
UV asymptotic solution towards small x, for a given input boundary condition T(0) the system has a 'preferred' 
dynamical curve corresponding to the powerlaw case T(0) = until 'forced' to deviate to the constant boundary value 

(ifr(o)^o). 

In the case of the pure powerlaw solution (numerically here, the lowest value of T(0) = 10 5 ), we find that the 
ghost propagator dressing function solution can be well reproduced over the entire momentum range by the form 



Dcc(x) 



1 



'(1 -a) 



(4.2) 



% /2X/VcY x ln(l + a;)' 

with a single parameter, a sa 0.689 (we find that a is slightly dependent on A, which has here the value A = 10 6 
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FIG. 3: Plot of the ghost renormalization coefficient Z c as a function of T(0). See text for details. 
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FIG. 4: Plot of the ghost propagator dressing function in the case T(0) = 
scales shown) compared with the fit formula Eq. (|4.2[) . See text for details. 



10 5 (indistinguishable from the powerlaw for the 



in the scaled units). The comparison is shown in Fig. Ul One can see that the form given in Eq. (J4.2I) will exactly 
reproduce the known infrared behavior because of the small x expansion of the logarithm. That the large momentum 
behavior is, up to a constant, given by \j\J\s\x was also found in Ref. [2J]. Recall that it was earlier seen that 
the propagator dressing function behaves as ~ x 7s for asymptotically large x: the 1/ylax behavior represents an 
alternative description to the naive perturbative anomalous dimension and indeed appears superior since it is valid 
over a wide range of momenta. 



A. Comparison with lattice results 

We have seen that using a lattice inspired result 34] as input for the spatial gluon propagator, the ghost propagator 
dressing function either has an infrared divergence characterized by the exponent 1/2, or is constant in the infrared 
depending on the boundary condition T(0). This result could be compared to the corresponding lattice results. 
Various studies |3ll433l |35| have reported an infrared divergence, i.e., a propagator dressing function of the form x~ K 
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FIG. 5: Plot of the ghost propagator dressing function in the case T(0) = 10~ 4 (indistinguishable from T(0) = on the 
scales shown) and T(0) = 0.2 compared with the powerlaws characterized by the exponents k — —j g (perturbative) , n = 0.2 
(illustrative lattice result) and k, — 1/2 (infrared analysis). The scale x = 0.25 roughly corresponds to the lowest lattice 
momenta. See text for details. 

for small x, but with an exponent k ss 0.2 — 0.25 (all studies agree on the Gribov form for th e eq ualtime spatial gluon 
propagator), apparently in contradiction to the results here and from the canonical approach [24|,[25| (where n = 1/2). 

The aforementioned lattice studies naturally involve a smallest infrared spatial momentum scale k 2 /m 2 = x m in 
and Xmin ~ 0.25 — 0.5. To extract the infrared exponent k, one requires a range of momenta for which the powerlaw 
relationship holds (or the constant T(0) is clearly visible). However, taking this range to be, for example, one order 
of magnitude (i.e., x = [0.25, 2.5]) and comparing with the results displayed in Fig. [21 one sees that all curves (with 
< r(0) < 0.3) are still within the transition region where the exponent changes from the perturbative value (x la ) 
to the infrared (a; -1 ' 2 ). In other words, it may be that the lattice results have not yet genuinely probed far enough 
into the infrared to see the full exponent. 

There is however the possibility that the lattice results might be infrared constant. In Fig. [5] as an illustrative 
example, we compare the results for the infrared divergent ghost propagator dressing function (numerical and infrared 
powerlaw, corresponding to k = 0.5) with the numerical curve for T(0) = 0.2 and a fit form with k = 0.2. Also shown is 
the perturbative asymptotic result. Importantly, for a wide range of values (x « [0.25, 10]) the curves with T(0) = 0.2 
and k = 0.2 are indistinguishable. Whilst this is only an illustration, it does give an alternative explanation to 
the lower lattice value of k as corresponding to a nonzero value of T(0) whilst not contradicting the earlier infrared 
analysis. Again, the deciding factor is that x m i n cannot be arbitrarily lowered. Note also that since the lattice results 
with x m in ~ 0.25 do see an infrared enhancement with k = 0.2 rather than a constant behavior for the dressing 
function, this would presumably place an upper bound on the value of the boundary condition: roughly T(0) < 0.2. 

To summarize, it would appear that in order to discriminate between the various infrared solutions, the lattice 
studies have not yet gone far enough into the infrared to give a clear signal that one has a particular value of the 
exponent k, or a particular boundary condition T(0). Such studies are planned [49( and using the analysis of the 
ghost Dyson-Schwinger equation presented in this study may be of help in the extraction of the infrared behavior. 



DISCUSSION OF THE RESULTS 



In this Section, we shall present and discuss arguments for the possible interpretation of the numerical results. To 
this purpose, we shall draw on various different ideas and concepts. In order to aid clarity, we split the discussion 
into different parts. The first part reiterates the identification of T(0) as a boundary condition. The second then 
addresses the issue of the physical interpretation of such a boundary condition and its connection to the Gribov 
gauge fixing ambiguity. The third part then continues by turning the arguments round, resulting in a conjecture that 
physical results may actually be independent of the value of the boundary conditions. The final part then discusses 
an alternate interpretation of T(0) as being concerned with the renormalization scheme and why this is seemingly 
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ruled out by the results here. 

1. r(0) as a boundary condition. 
That the various infrared solutions of the ghost Dyson-Schwinger equation are delineated by some constant, r(0) 



should not come as a surprise. Indeed, this has been explicitly seen in the canonical approach to Coulomb gauge 27 1 
and there is a similar situation in Landau gauge as discussed in the Introduction. The Dyson-Schwinger equations 
are derived as functional differential equations that can only relate Green's functions and as such, require (functional) 
boundary conditions to finally specify a value for a given Green's function. Clearly the solutions should agree with 
the asymptotic perturbative expansion in the vicinity of the renormalization point and the functions should be 
smooth and continuous for finite, non-zero momenta (a feature naturally present for an iterative numerical solution). 
The remaining boundary condition in the case of the ghost Dyson-Schwinger equation is then the value of T(0). 
Normally, one would be considering a coupled system of Dyson-Schwinger equations whereby the solution for the 
gluon propagator is also being sought and the issue of the existence of the boundary condition is invariably obscured 
by technical details (such as the truncation scheme); here we have utilized the general consensus about the Gribov 
form for the equaltime spatial gluon propagator to fix this sector so that the existence of the boundary condition as 
input to choosing the relevant solution is made explicit. Indeed, this is particularly relevant in Coulomb gauge within 
the functional integral formalism, where the energy divergence of the ghost-loop and the necessary nonperturbative 
cancellation prohibits a direct back-coupling of the ghost to the gluon polarization. 

2. T(O) and the Gribov problem. 

Having established the requirement for a boundary condition, the important question is thus: is there a particular 
physical value? To answer this question we obviously need to know to what physics the different solutions correspond 
to. The most pertinent aspect of the numerical results presented here is that for a wide range of values (0 < T(0) < 0.3), 
the UV solution and the ghost renormalization coefficient are unaltered and meaning that the boundary condition 
has nothing to do with either the perturbative region or the renormalization condition (at least conventionally where 
the renormalization scale is chosen to be perturbative). We can thus conclude that the boundary condition, at least 
for the ghost Dyson-Schwinger equation in isolation, is connected only to the nonperturbative (infrared) physics. 
Additionally, that the various solutions exhibit a particular infrared behavior where as one goes from large to small 
momenta, the curves follow the dynamical solution until 'forced' to deviate to the constant boundary value suggests 
that the boundary condition is connected to a freezing of the dynamics below a particular scale. 

Let us now introduce some aspects of Coulomb gauge and the Gribov problem relevant to the discussion here. In 
the functional approach to Coulomb gauge (within the first order formalism, but the same applies here) , the system 
can be shown to reduce to the 'would-be-physical' degrees of freedom QQ, [39|, |40|: that is, after resolving Gauss' 
law only two transverse spatial gluon degrees of freedom remain (in quantum electrodynamics these would be the 
physical photon polarization states) and the longitudinal, temporal, and ghost degrees of freedom cancel. In the 
Gribov-Zwanziger confinement picture |14l4l6l |. the transverse spatial propagator (given by the Gribov formula here) 
is infrared suppressed and subsequently drops out of the physical spectrum. The reduction to transverse spatial 
degrees of freedom has the Faddeev-Popov (FP) operator as a central element and the ghost propagator is precisely 
the expectation value of the inverse of this operator. However, the known existence of zero modes of the FP operator 
and the associated Gribov ambiguities complicate matters [14J . The Gribov problem is a result of the fact that 
the gauge is not completely specified. The lack of complete gauge fixing manifests itself in gauge- variant quantities 
(i.e., the propagators) and implies an ambiguity in the definition of the underlying functional integral. Given that 
the spatial gluon propagator in Coulomb gauge corresponds to the 'would-be-physical' degrees of freedom (in the 
above sense), one might expect that the Gribov problem would be relatively unimportant - we have, after all, the 
Gribov form as input; on the other hand, the ghost propagator would be entirely dependent on the resolution of the 
incomplete gauge- fixing (explicitly found to be the case in 1+1 dimensions [50]). This scenario is the case here with a 
fixed spatial gluon propagator (the Gribov form) and a ghost that must be specified with further input but that does 
not directly back-couple to the gluon. 

The first part of the argument for a connection between the Gribov problem and the boundary condition T(0) goes 
as follows. In the presence of the Gribov ambiguity, the form of the Dyson-Schwinger equations (at least in more 
than one spatial dimension for Coulomb gauge) is unchanged [51| . The implication is that the dynamical behavior 
of the solutions would not be affected by the existence of Gribov copies and this is consistent with the characteristic 
pattern of the results here (Fig. [5]) where with various values of T(0) (at least for the range < T(0) < 0.3) as T(0) is 
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lowered, the solutions follow the same dynamical curve further into the infrared before freezing out to the value T(0) 
imposed as the boundary condition. 

The second part of the argument connecting L(0) and the Gribov ambiguity concerns the following observation. 
Generally, when (improperly) integrating over the gauge group, gauge-dependent functional integrals (propagators) 
would be suppressed by averaging over the gauge copies. This is an important point: gauge copies either have 
the same phase in group space when contributing to the functional integral in which case the integration over the 
gauge copies is merely an unobservable constant in the normalization; or if the gauge copies have a different phase 
then the integration over the group space can only reduce a gauge-dependent functional integral. Gribov copies are 
accompanied by zero modes of the Faddeev-Popov operator, which for Coulomb gauge reads 

-V-D ab [A], D ab [A}=5 ab \7-gf acb A c {x). (5.1) 

Perturbatively, in the sense of an expansion in the gauge field (A) around zero there are no Gribov copies such that the 
perturbative regime would be insensitive. The ghost propagator results are insensitive to T(0) at large momenta where 
perturbation theory is valid (as explicitly demonstrated previously). In the nonperturbative case, for low momenta, 
there are 'more' Gribov copies (we do not know how to count the number, but from the form of the Faddeev-Popov 
operator, the large fields A within the covariant derivative that account for the existence of Gribov copies are large 
in the sense of their relationship to the spatial differential operator, which in momentum space translates, loosely 
speaking, to the spatial momentum) . There are two plausible signals for the possible presence of Gribov copies within 
the functional integral for the ghost propagator arising from improper gauge averaging: 1) below a certain momentum 
scale, the averaging over the gauge group induced by the existence of Gribov copies might serve to freeze the functional 
integral to a specific value and 2) for successively lower momenta, the difference between a functional integral affected 
by Gribov copies and one that isn't would increase. The results in this study show that for finite r(0), there is 
precisely such a freezing below a particular scale whereas for T(0) ~ there isn't. Also, the difference between the 
various solutions with finite and zero values of T(0) does increase with successively lower momenta. Now, T(0) is a 
boundary condition specifying a particular allowed solution of the Dyson- Schwinger equation and thereby implicitly 
selecting a definition of the ghost propagator functional integral from the various possibilities given that there is an 
ambiguity in the gauge fixing. That for different values of r(0), the results explicitly exhibit the characteristics that 
one would expect in the possible presence of Gribov copies suggests that T(0) is indeed intimately connected to the 
Gribov problem. The case r(0) = would appear to correspond to that implicit definition of the functional integral 
for the ghost where the Gribov copies have no effect, i.e., the curve is purely dynamical and there is no freezing in the 
infrared. For finite r(0), the dynamical curve freezes below a certain scale indicating the influence of Gribov copies 
and their averaging. 

There are two further arguments connecting T(0) and the Gribov problem that are known in Coulomb gauge: 
Gribov's original work [lj] and the case of 1 + 1-dimensions [50] ■ In Ref. [lj] , the perturbative aspects of the ghost 
propagator were discussed. In the language used here, in order to ensure the absence of gauge copies (in the sense 
that the FP operator be positive definite), the restriction T(0) > was introduced. In the case of 1 + 1-dimensional 
Coulomb gauge (where because the 'physical' transverse spatial propagator cannot exist, the system reduces to purely 
a discussion of the gauge-fixing and Gribov copies) it was seen explicitly that the ghost propagator functional integral 
after resolving the Gribov problem was that with effectively the lowest allowed positive value of the boundary condition 
(due to the spatial compactification of the underlying manifold employed there, one could not assign a value of T(x) 
for x = 0). Functional integrals constructed to include Gribov copies had different boundary conditions, but obeyed 
the same Dyson-Schwinger equation. 

One objection to the above assertion that T(0) = corresponds to the situation where Gribov copies have no 
effect is of course that the lattice data may be consistent with the solution for T(0) ~ 0.2 as discussed earlier. That 
the ghost dressing function is infrared finite seems to be the case in Landau gauge on the lattice (see for example 
Ref. [l0| and references therein). However, it could be that the different results are not in conflict but are results 
for different problems. Gribov copies are the consequence of an incomplete gauge fixing, which manifests itself in 
the eigenspectrum of the Faddeev-Popov operator and this eigenspectrum is obviously dependent on the boundary 
conditions one assigns to the eigenfunctions. Namely, the Gribov problem may be different for a compact manifold 
(the lattice with periodic boundary conditions) or in the continuum (with a flat Minkowski/Euclidean metric and 
boundary conditions at infinity) and also different when comparing Coulomb and Landau gauges [52| . When inserted 
into the functional integral, it is then currently an open question as to what effect this might have. 

3. r(0) as an implicit gauge choice. 

Another way to view the connection between T(0) and the Gribov problem is that it can be conjectured that 
the additional (nonlocal) gauge fixing required to eliminate Gribov copies is linked to the discussion of a specific 
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choice for T(0). However, turning the argument around, one could interpret T(0) as a choice of gauge, as was 
proposed in Ref. [5J|. More properly, the interpretation of T(0) would be as a choice of gauge completion: the overall 
gauge choice (Landau or Coulomb gauge, for example) is required to have a well-defined perturbative propagator 
and because the dynamical solution to the Dyson-Schwinger equation must be smoothly and continuously connected 
to this perturbative expression in the ultraviolet, then the lack of complete gauge fixing does not mean that the 
Green's function is ill-defined, rather that there exist multiple solutions and the choice of gauge completion defines 
the nonperturbative component of the propagator in the sense that a specific solution (corresponding to a particular 
value of the boundary condition) is chosen. Indeed, the interpretation of the boundary condition T(0) as an implicit 
choice of gauge is rather convenient since as long as one can show that the different choices are physically equivalent, 
one can choose the gauge so as to make specific calculations more straightforward. That is to say that whilst the 
elimination of Gribov copies may indeed provide a definite value for T(0), the arbitrariness is to do with the gauge 
fixing and one may specify any gauge one desires. This provides a possible answer to the original question about the 
existence of a particular physical value for the boundary condition T(0) - there may be no unique physical value, just 
as propagators are different in different gauges. Furthermore, issues such as the confinement mechanism may appear 
different for different choices of T(0). The task would of course be to eventually show that the same physical results 
can be obtained with different boundary conditions. We shall briefly discuss the case T(0) = below. 

4- r(0) as a renormalization scheme. 

As a final remark here, another interpretation of the boundary condition is that the value of T(0) is a choice of the 
renormalization scheme to absorb constants into the renormalization coefficient Z Cl and which again is a free choice 
(and which has no impact on the perturbative renormalization). The Gribov-Zwanziger confinement scenario [14j - 4l6| 
motivates the choice T(0) = 0. However, that the renormalization coefficient Z c is seen to be stable for a range of 
values, T(O) (see Fig. [3]), would seem to rule out this interpretation of Z c for this aspect of the Gribov-Zwanziger 
scenario in Coulomb gauge. 

VI. CRITICAL GHOST CONFINEMENT AND INFRARED SLAVERY 

As has been seen, when the boundary condition T(0) = 0, the ghost propagator dressing function behaves like l/\k\ 
in the infrared. This is of course the result for the so-called critical solutions obtained in the canonical formalism 
and provides a very simple picture for many aspects of confinement 28, _29| . In the functional formalism, the infrared 
critical (powerlaw) ghost solution provides an important possible addition to this list, as we shall now motivate. 

The ghost and spatial gluon propagators have been discussed so far, but in the functional approach to Coulomb 
gauge there is a third propagator - the temporal gluon propagator corresponding to the temporal component of 
the gluon field (a = A ). In the canonical approach to Coulomb gauge, one starts with Weyl gauge (tr = 0) and 
then imposes Coulomb gauge [2J, [30] , so the temporal gluon propagator does not exist in this sense (the analogous 
quantity to the temporal gluon propagator in the canonical formalism is the non-Abelian color potential) . Within the 
functional formalism in the heavy quark limit, it has been shown that when the temporal gluon propagator diverges 
like l/k (or equivalently, that the dressing function goes like 1/fc 2 ) in the infrared, one directly obtains a linear rising 
potential (and bound states) for color singlet quark-antiquark and diquark (for N c = 2 colors) configurations with no 
finite energy colored states permitted and corresponding precisely to the old infrared slavery picture of confinement 
[44| . The lattice results for the temporal gluon propagator do exhibit this infrared behavior I3lll3a | ■ 

The decomposition of the temporal gluon propagator and proper two-point function reads [4l| : 

w£(k)=i6 ab Daai f°' k2) , Tf a (k) = -zS ab k 2 T (7<T (k 2 ,k 2 ), D a(J (kl 1 P)Y aa {klP) = l. (6.1) 

k z 

In addition, the decomposition of the proper spatial gluon two-point function (the spatial gluon polarization) is 



rf A a(k) = -is ab 



k%(k)r AA (k 2 ,k 2 ) + k 2 ^r AA (k 2 ,k 2 ) 

k 2 



(6.2) 



where the dressing function T AA is the inverse of D AA and T AA is the longitudinal dressing function component of the 
polarization (in Coulomb gauge, this quantity is emphatically not vanishing unlike in Landau gauge). Importantly, 
the Slavnov- Taylor identity for the two-point functions tells us that [431 ] 

T aa (kl k 2 ) = T AA (kl k 2 ) \T- cc {k 2 )] 2 (6.3) 
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(we replace the subscripts for Tc C to avoid confusion). On general grounds, it can be shown that D aa (and consequently, 
Tgn-) must have a nontrivial instantaneous (energy independent, or ~ &(x\ — x%) in configuration space) component 
48]. In fact, lattice results indicate that there is only a very weak energy dependence for D aa [35| and so it appears 
that neglecting the non-instantaneous component is a good approximation. In other words 

Ara(*§, k 2 ) = D aa (k 2 ), T aa {kl k 2 ) = T aa {k 2 ), T AA (k 2 , k 2 ) = T AA (k 2 ). (6.4) 

Assuming multiplicative renormalizability and assigning the renormalization coefficient \J Z a to the temporal gluon 
field a as before, since both the longitudinal and transverse components of the spatial polarization share the same 
renormalization coefficient, then we immediately see that 

Z a = Z A Z 2 (6.5) 

as a consequence of the Slavnov- Taylor identity, as has been known for some time [lf|. With Eq. (|2.13[) . the product 
g 2 D aa is thus a renormalization group invariant, just as g 2 D AA D^ c is. Taking the infrared powerlaw form for the 
ghost propagator dressing function (noting the scaling of the variables with to 2 ), one infers that 

g 2 D <Ta {k 2 ) P ^° ^_ m \ (6.6) 

k 2 T AA {k 2 ) 

meaning that if T AA is a constant in the infrared, one has exactly the temporal gluon interaction required for a linear 
rising potential in Coulomb gauge. Recall that the lattice results for the spatial gluon propagator indicate that the 
Gribov scale to is not dependent on the renormalization scale [34] and one had the approximate relationship to ~ 2^/a w 
where u w is the Wilsonian string tension. The numerator factor to 2 arising from the Gribov formula would be thus 



entirely consistent with the coefficient of the linear rising potential discussed in Ref. [44J . Also recall that in Coulomb 
gauge, the spatial gluon propagator corresponds to the 'would-be-physical' transverse degrees of freedom which indeed 
suggests that the Gribov scale should be related to a physically observable scale. So, if T AA were to be a constant in 
the infrared, there would be a direct connection between the spatial gluon propagator and the confinement potential 
even though the spatial gluon propagator is infrared vanishing (as suggested by the Gribov-Zwanziger confinement 

scenariojfj-G]!)- 

That T AA is a constant can be argued by appealing to the dimensional type of infrared analysis employed earlier 
for the ghost Dyson-Schwinger equation and general arguments for Dyson-Schwinger results seen in Landau gauge. 
In terms of renormalized quantities and recognizing that the renormalization coefficients are dependent on the regu- 
larization scale, A, the ghost Dyson-Schwinger equation reads schematically 

Tg(k 2 ;^g) = Z c ( f i,g;A) + Z A (^g ] A)^- c (k 2 ;^g;A). (6.7) 

Now, with the result that Z A = 1, the self-energy integral must separate into two components, i.e., 

^- cc (k 2 ;fi,g;A)^^(k 2 ; l i,g) + ^( f i 1 g;A) (6.8) 

such that the equation can be renormalizable. The earlier dimensional analysis for the infrared could then be performed 
without reference to the (large) scale A and S|*" after a straightforward subtraction of the equation for k 2 = \i. Then, 
counting the powers of k 2 led directly to the powerlaw exponent if one assumed the absence of the infrared constant 
solution and disregarded the Gribov scale to. In the case of the Dyson-Schwinger equation for T AA , one has the 
similar schematic form 

T AA {k 2 ; n,g) = Z A ([i,g; A) + £ Z 2 ( M , g; A)£,(fc 2 ; /*, g; A) (6.9) 

i 

but with many more loop self-energy terms (labeled with the index i). The renormalization coefficients Zi are the 
coefficients for the dressed vertices that enter the corresponding loop integrals and can be explicitly derived from 
the Slavnov- Taylor identities j43|. The important information is that only the spatial ghost-gluon vertex coefficient 
is trivial (Z A = 1). Thus, only the ghost loop of the Dyson-Schwinger equation for T AA has a form similar to the 
ghost self-energy above, as is also the case in Landau gauge. However, in Coulomb gauge this loop is explicitly 
energy divergent and must cancel. Therefore, one cannot separate the finite (A-independent) and divergent (fc 2 - 
independent) parts of any of the loop integrals Sj and the dimensional analysis, which relies on this decomposition, 
cannot hold. This argument would also be true for the Dyson-Schwinger equations for T AA (the transverse part of 
the polarization). Actually, the lack of such a separation between infrared and ultraviolet parts of the self-energy 
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terms is another of the difficult technical problems associated with the functional approach to Coulomb gauge and is 
an explicit demonstration of the difference to Landau gauge (where it is the ghost loop that is the important infrared 
contribution to the polarization Q). 

In the absence of the dimensional analysis, the infrared character of the Dyson-Schwinger equations cannot be 
directly inferred. However, one can use the experience gained from Landau gauge to anticipate the likely outcome. It 
has been known for a long time that in Landau gauge, there exist potential quadratic UV divergences in the tensor 
component of the polarization that is proportional to the metric, which here would be the transverse part of the 
spatial gluon polarization (i.e., the Dyson-Schwinger equation for Taa) [5J| (see also Ref. [22} and references therein). 
Thus, one would naively expect terms in the unrenormalized equation for the dimensionless dressing function that 
behave like A//c 2 , but which should cancel with a proper regularization procedure (and ignoring in this case the energy 
dependence). After regularization and renormalization (which doesn't alter the /c 2 -dependence), such terms can only 
be of the form m 2 /k or ji/k 2 , with any other terms being infrared subleading. Thus, one would expect an infrared 
behavior for Daa — ^aa °f ^ ne f° rm 

Daa ^ ° k 2 , (6.10) 

which is exactly the case for the Gribov formula (and observed on the lattice). For the longitudinal part of the 
polarization, one would expect logarithmic UV divergences and in the presence of the additional Gribov scale m, this 
would indeed favor a constant type of infrared solution, just as is the case for the dressing functions in the fermion 
sector of QCD (see for example Ref. [55] ) : 

Taa ~ const. (6-11) 

In other words, one can justify that Taa is indeed constant in the infrared such that with the ghost boundary condition 
F(0) = 0, the temporal propagator would be precisely that required for a linearly rising confinement potential. 

Obviously, the discussion of this section is somewhat speculative and is certainly intended more to discuss the 
plausibility of a simple confinement scenario in Coulomb gauge than as a quantitative calculation. One example of a 
technical complication that might arise in such a quantitative study is that the Slavnov- Taylor identities are in general 
deformed by the presence of a UV cutoff scale A: in Landau gauge, this has been explicitly studied in Ref. J22J (and 
references therein). However, what should be amply clear is that the confinement mechanism in Coulomb gauge is 
certainly dependent on the ghost boundary condition F(0), underscoring the necessity of pursuing the understanding 
of the ghost propagator in the infrared further. 

VII. SUMMARY AND OUTLOOK 

In this study, the ghost Dyson-Schwinger equation in Coulomb gauge has been considered numerically using lattice 
input for the spatial gluon propagator. It is demonstrated that the dynamical solution must be explicitly supplemented 
by a nonperturbative boundary condition (the value of the inverse ghost propagator dressing function at zero spatial 
momentum). With various values of this boundary condition, the solutions have a characteristic behavior: in the 
ultraviolet, all solutions lie on top of each other and agree with the asymptotic leading order perturbative result; going 
from high to low momenta, the solutions follow the common dynamical curve into the infrared until freezing out at 
a point determined by the boundary condition. In this way, it is seen that both critical (powerlaw) and subcritical 
(infrared finite) solutions exist. It was seen that the critical solution is characterized in the infrared by the exponent 
a = 1/2. The ghost renormalization coefficient was observed to be largely independent of the boundary condition. A 
qualitative comparison to lattice results for the ghost propagator dressing function was illustrated, with the conclusion 
that it would be desirable to have lattices that go further into the infrared in order to discriminate between the critical 
and subcritical solutions. The characteristic pattern of the solutions for different boundary conditions can be naturally 
interpreted in terms of the Gribov gauge fixing ambiguity and it is conjectured that if the boundary condition is related 
to a choice of gauge fixing completion, then physical (gauge invariant) results should not be dependent on the choice. 
Plausibility arguments are put forward that connect the critical solution for the ghost propagator dressing function 
with the temporal gluon propagator and the confinement mechanism. 

Whilst the aim of this study was to demonstrate the explicit requirement for a nonperturbative boundary condition 
to supplement the dynamical solution of the ghost Dyson-Schwinger equation in Coulomb gauge and to present 
numerical results, the characteristic behavior that was found led to two primary conjectures. Clearly then, the task 
will be to investigate further. The first conjecture is that the boundary condition studied here is connected to the 
resolution of the Gribov problem in Coulomb gauge. Work is currently underway to study this connection explicitly 
[56j. The second conjecture concerns the connection of the ghost propagator to confinement via the temporal and 
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longitudinal spatial gluon proper two-point functions courtesy of the Slavnov- Taylor identities. This is a formidable 
task because of the inherent technical difficulties of solving the Dyson-Schwinger equations within the noncovariant 
Coulomb gauge functional formalism; however, these technical problems are being steadily overcome and one may 
hope for concrete results within a reasonable timeframe. 
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